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We show the appearance of spatiotemporal stochastic resonance in the Swift-Hohenberg equation. 
This phenomenon emerges when a control parameter varies periodically in time around the bifur- 
cation point. By using general scaling arguments and by taking into account the common features 
occurring in a bifurcation, we outline possible manifestations of the phenomenon in other pattern- 
forming systems. 



PACS numbers: 05.40. -fj, 47.54-fr 

In the last years the phenomenon of stochastic reso- 
nance (SR) l|-||] has been the subject of intense activ- 
ity, to the extent of becoming cross-disciphnary due to 
the great number of apphcations to different fields of sci- 
ence. The main result of SR, which in some ways could 
be considered counterintuitive, shows a constructive role 
of noise, since the response of a system to a periodic sig- 
nal may be enhanced with the addition of an optimized 
amount of noise. In this sense, the presence of SR is well 
characterized by the appearance of a maximum in the 
output signal-to-noise ratio (SNR) at a certain nonzero 
noise level. 

In spite of the fact that spatial heterogeneity is one of 
the most common features of systems away from equi- 
librium, up to now there has not been a complete un- 
derstanding about the phenomenon of SR in spatially 
extended systems. In this context, only a few recent 
results are available from the literature M-Q. A re- 
markable feature of spatial systems is their possibility 
to develop patterns that can form via bifurcations, for 
instance from the spatially uniform state, as a control 
parameter is varied p^,0. In this situation, stochas- 
tic perturbations play a crucial role in the initial stages 
of pattern formation to the extent that the system may 
display macroscopic manifestations of thermal noise. 

In this Letter we will show that temporal periodic vari- 
ations of the control parameter in the presence of noise 
can lead to spatiotemporal stochastic resonance (STSR) 
when the system is close to a bifurcation point. To be 
explicit, we will treat the Swift-Hohenberg (SH) equa- 
tion ||l^,|l6|,0 in detail, which models Rayleigh-Benard 
convection near the convective instability. This exam- 
ple, which constitutes the paradigm of pattern forma- 
tion, exhibits many features common to systems around 
a bifurcation point. In this regard, and for the sake of 
generality we will outline applications of our results to 
other situations. 

In the vicinity of the instability point, Rayleigh-Benard 
convection can be described by means of the stochastic 
SH equation, which in dimensionless spatial units is given 
by 
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Here the control parameter h{t) = — k + asn\{LOQt), with 
K, a and wq constants, reflects the presence of an ex- 
ternal periodic forcing, due, for instance, to variations 
of the temperature difference between the plates of the 
convective cell [ |lO|Jl7|Jlq] . Moreover, g and q are param- 
eters which depend on the characteristics of the system 
and ^(r, t) is Gaussian white noise with zero mean and 
second moment {£_{r,t)£,{r',t')) = 2D5{r - f')5{t - t'), 
defining the noise level D. 

In order to describe the temporal evolution of the sys- 
tem, we will consider the convective heat flux, which in 
this model is given by 



J{t) =c / il:{r,tfdr 



(2) 



where c is a constant depending on the physical charac- 
teristics of the system. This magnitude constitutes the 
order parameter of the transition from the homogeneous 
state to the state where spatial structures develop. In 
regards to spatial order, it can be revealed by the time- 
averaged structure factor 



S{k) = lipki'-k 



(3) 



where ipk is the spatial Fourier transform of the field tp 
and < >t indicates time and noise average. A sharp 
peak in this magnitude makes the presence of an ordered 
spatial structure manifest. 

Since we are interested in the effects of noise, we will 
first analyze how a small amount of noise affects our sys- 
tem. In the absence of noise, irrespective of the initial 
condition, the field ip goes to zero at large times. For 
a sufficiently low noise level, and far from the possible 
initial transient, ip is small and the nonlinearity in Eq. 
(nf) does not play any role. In this situation, by using 
dimensional analysis it is easy to see how the character- 
istic magnitudes scale with noise. We note that, when 
the linearized equation is considered, any dimensionless 
parameter cannot depend on the noise level, because only 
D involves the dimensions of the field tp. Thus tp scales 



with the noise a.s 'ip (x \JD. The SNR has dimensions of 
the inverse of time B, then for J(t) it is given by 



SNR = c^o/i (7) 



(4) 



where /i is a dimensionless function which depends on 
the set of dimensionless parameters h/ojq, a/ujQ and q/loq, 
denoted by 7. We then conclude that the SNR does not 
depend on the noise level. However, both the signal and 
output noise scale with the noise as D^. This fact in- 
dicates that the output signal increases when noise in- 
creases. Additionally, the structure factor also follows a 
scaling law: S{k) oc D. 

Since for low noise level the SNR does not depend on 
D, the lowest order correction in D to the constant value 
of the SNR comes from the nonlinear term. To eluci- 
date its form, it is convenient to rewrite Eq. (||) in the 
following way 



'di 



-J- = [-{k + 5V') + a sin(woO - 9(1 + V')']^' + ^(r, t) 



(5) 



Due to the fact that for low noise level tp'^ oc D, in first 
approximation k + gip'^ can be interpreted as an effective 
parameter k = k + gujQ f2{"f)D, with /2 being a positive 
dimensionless function. Consequently, Eq. dq) has the 
same form as the linearized version for which the scaling 
law for the SNR [Eq. fl)] has been derived. By replacing 
K by K in Eq. (0) and expanding around k we then obtain 



SNR « uofi (7) + 1^5^0^72(7)^ 

OK 



(6) 



which includes the lowest order correction in D to the 
SNR due to the nonlinear term. An important conse- 
quence of this result is the fact that the knowledge of 
the dependence of the SNR on the parameter k, for the 
linearized equation, enables us to predict the presence of 
SR when the nonlinear term is considered. If /i is an 
increasing function of k then dfi/dn is positive and con- 
sequently the SNR is an increasing function of D for low 
noise level. Since for high D the SNR decreases, one then 
concludes that it exhibits a maximum thus indicating the 
presence of SR. 

To verify these results we have integrated the previ- 
ous equations by discretizing them on a mesh p9| and 
then by using a standard method for stochastic differen- 
tial equations [g^ . For the sake of simplicity, we will first 
consider the one-dimensional case, although the previous 
scaling laws hold independently on the dimensionality of 
the system. In Fig. |^(a) we have depicted the SNR as a 
function of k for particular values of the remaining pa- 
rameters. One can see that this quantity has a maximum 
at Kjnax- As a consequence, in view of Eq. (^) the SNR 
may increase (k < Umax) or decrease (k > Kmax) with 
D. In Figs. 0(b) and |l](c) we have represented the SNR 



as a function of D for n < Kmax and k > Kmax, respec- 
tively. These figures corroborate the dependence of the 
SNR on the noise level and on the parameter g in Eq. 
(^). This scaling argument is robust upon varying the 
nonlinear term. For example, if one replaces the term 
gtp"^ by gIV'IV' one obtains a SNR that increases as a/D 
instead of as D. This result is shown in Fig^(d). In this 
regard, other nonlinearities have also been successfully 
tested. 

The previous analysis has shown the possibility for the 
appearance of SR in the SH equation. In Fig. ||(a) we 
have plotted the SNR for J{t) as a function of D, ob- 
serving a maximum at a nonzero noise level. It is worth 
pointing out that an optimized amount of noise increases 
the SNR up to 20 dB. As far as the spatial effects are 
concerned, the sharpness of the peak of S{k) can be an- 
alyzed by considering its height over its variance, which 
will be denoted by Q. This magnitude has also a max- 
imum [Fig. |2|(b)] which is due to the fact that for low 
noise level Q scales with D in the same way as S{k), 
whereas for high noise level the noise destroys the spatial 
structure. An example of the spatiotemporal evolution 
of the system is depicted in Fig. ||. For sufficiently low 
noise level, the system exhibits neither spatial nor tem- 
poral structures. For intermediate values of the noise, 
the system shows a spatial pattern and a coherent re- 
sponse to the periodic variations of h{t). These patterns 
are destroyed for higher values of D. In this case, noise 
induces periodic spatial patterns and ordered temporal 
behavior, thus playing a constructive role in both space 
and time. 

The characteristics observed for the one-dimensional 
SH equation also hold in two dimensions. To illustrate 
the effect of noise in the two-dimensional case, we have 
plotted three patterns for different values of D in Fig. |j. 
It becomes clear that there exists an optimum noise level 
in order to observe the typical convective rolls. 

Our results can also be applied to a great variety of 
systems in the vicinity of a bifurcation point. We note 
in this context that the form of the scaling laws does not 
depend on the way in which the nabla operator enters 
the linearized equation. A remarkable example, which 
enables one to predict the presence of STSR by only con- 
sidering scaling arguments, is the generalized Ginzburg- 
Landau equation [gl],|0|. In addition to the cubic term 
discussed previously, this equation presents a quintic one 
responsible for saturation effects. Since in this situation 
the coefficient of the cubic term may be positive or nega- 
tive, there will be a case in which the SNR is an increasing 
function of the noise level giving rise to the appearance 
of STSR. 

In summary, for the first time we have shown the pres- 
ence of STSR arising when the system undergoes a bifur- 
cation, a usual mechanism for pattern formation. One 
aspect to be emphasized is the fact that in spatial non- 
linear systems, the presence of noise in combination with 



a periodic signal may give rise to ordered spatiotempo- 
ral structures which are not present in the absence of 
noise. We have also explicitly shown the presence of 
STSR in the SH equation, due to its great importance 
as a model of pattern formation. However, our results 
can be straightforwardly applied to a wide variety of sys- 
tems, provided they exhibit the common characteristics 
of this bifurcation. These findings therefore open up new 
perspectives about the general consideration of the phe- 
nomenon of SR to the field of pattern-forming systems. 
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FIG. 1. (a) SNR from the linearized SH equation as a 
function of k for a = 1, g = 1, ujo/2-k = 0.195, and g = 1. (b) 
ASNR = SNR - SNR(i3 = 0) from Eq. (Q) as a function of 
D for K = 0.3, g = 1, ujo/2-k — 0.195, and g = 0.5 (circles), 
<7 = 1 (squares) and g — 2 (triangles). The points have been 
fitted by a power law, obtaining an exponent equal to 0.99 for 
the three values of g. (c) ASNR as in case b) but for k — 3. 
In this case the exponents of the power law are 0.99, 1.00, 
1.02. (d) ASNR from Eq. mh when replacing the nonlinear 
term by gVI'^l) ^■s a function of D for k = 0.3, q — 1, and 
ujo/2tt = 0.195, and g = 1. In this case the scaling exponent 
is 0.50. In all figures the size of the system is 32. 
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FIG. 2. (a) SNR and (b) Q as functions of the noise level 
D for a = 1, K = 0.1, ff = 1, g = 1, and cJo/27r = 0.024. The 
system size is 32. 



FIG. 3. Spatiotemporal evolution of the one-dimensional 
SH equation with periodic boundary conditions for the noise 
levels (a) D = 0.5 x 10-^^ (b) D = 0.5 x 10"*, (c) 
D = 0.5 X 10"^, and (d) D = 0.5. The values of the re- 
maining parameters are a — 1, k — 0.1, g — 1, q — 1, and 
ijj(,l2-K = 0.024. The system size is 32. We have also indi- 
cated the time evolution of the control parameter h{t). For 
each pattern, black and white colors stand for minimum and 
maximum values, respectively. 




FIG. 4. Representation of the field tp for the 

two-dimensional SH equation with periodic boundary con- 
ditions. The noise level is (a) D = 0.25 x 10~^°, (b) 
D = 0.25 X 10~^, and (c) D = 2.5. The remaining values 
of the parameters in all cases are a — 1, k — 0.1, g = 1, 
q — 1, and aJo/27r — 0.012. The system size is 50 x 50. 



